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This paper presents an integrated approach to the modeling and control of aircraft with 
flexible wings. The coupled aircraft rigid body dynamics with a high-order elastic wing 
model can be represented in a finite dimensional state-space form. Given a set of desired 
output covariance, a model reduction process is performed by using the weighted Modal 
Cost Analysis (MCA). A dynamic output feedback controller, which is designed based on 
the reduced-order model, is developed by utilizing output covariance constraint (OCC) 
algorithm, and the resulting OCC design weighting matrix is used for the next iteration 
of the weighted cost analysis. This controller is then validated for full-order evaluation 
model to ensure that the aircraft’s handling qualities are met and the fluttering motion of 
the wings suppressed. An iterative algorithm is developed in CONDUIT environment to 
realize the integration of model reduction and controller design. The proposed integrated 
approach is applied to NASA Generic Transport Model (GTM) for demonstration. 


I. Introduction 

New energy efficient aircraft concepts are being studied by NASA, industries, and academia, through 
better aerodynamic design, use of light weight structures, and novel flight control concepts. In particular, 
lightweight aircraft design has attracted considerable attentions in recent years in an effort to improve cruise 
efficiency and lower the induced drag. As structural flexibility increases, aeroelastic interactions with aircraft 
aerodynamic forces and moments can have adverse impact on aircraft’s stability and performance. More 
specifically, the modal frequencies of flexible wings can be in the range of aircraft’s rigid body modes and 
within the flight control bandwidth. Therefore, designing an effective aircraft flight control law becomes a 
challenging task. 

The aeroelastic aircraft control problem considered in this paper utilizes work of two disciplines; control of 
flexible structures and aircraft flight control law development. In the context of flexible structure controls, one 
must first develop an adequate mathematical model, which is reduced in size, prior to designing an suitable 
control. A number of model reduction methods 1-5 can be used for carrying out the order reduction. One 
method that we find it most appropriate for our study is the Modal Cost Analysis . 5 In MCA the total cost 
is expressed as a sum of modal cost of each aeroelastic mode to the defined output covariance cost function. 
Therefore, by examining each modal cost, we can attain a reduced-order model by retaining those modes 
with higher modal cost. The performance requirements for vibrational aeroelastic wings subject to random 
gust load can be prescribed in terms of output covariance. Since both model reduction and aeroelasticity 
effect use covariance as a measure of preformance, we utilize the output covariance constraint (OCC) control 
algorithm 9 to develop an integrated modeling and control algorithm to suppress the vibrational motion of the 
wings. The essence of OCC problem is finding appropriate weighting matrices so that a stabilizing controller 
is obtained while minimizing the control effort and meeting constraints on the output covariance. 

In the context of designing a flight control law, we utilize CONDUIT (CONtrol Designer’s Unified In- 
terface ), 11 a state-of-the-art multidisciplinary computational facility, in attaining optimal solution to the 
aeroelastic aircraft control problem. In CONDUIT, both handling qualities specifications for aircraft rigid 
body dynamics and OCC for aeroelastic wings are evaluated throughout optimization process. In this case, 
weighting matrices are considered design parameters to be tuned at every CONDUIT simulation. There- 
fore, the aeroelastic aircraft control problem, consisting of model reduction, OCC control, and flight control 
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design, can be entirely setup in CONDUIT environment. A preliminary analysis result for NASA Generic 
Transport Model (GTM) is presented. 


II. Coupled Aeroelastic and Aircraft Dynamics 

The coupled aircraft rigid body dynamics with 2 N aeroelastic wing modes; N bending modes and N 
torsional modes, can be described by 


Af e x e T C e x e T K e x e T C a x a [Dy Bs\ 


S s 


x a — A a x a T A s x e T .Af/Xg T B e S e T \Bj 


5 / 

<5 S 


( 1 ) 


where x e = [u>i W2 ■ • • wn #2 • • • € R 2N denotes the displacement for bending and torsional flexible 

modes, x a € R Ua the aircraft’s rigid body states, e.g. angle of attack and pitch rate for short-period mode, 5 e 
denotes the elevator/aileron/rudder deflection angle, <5/ and 5 S denote, respectively, flap and slat deflections. 
The size of x e can be very large, hence we consider the model given in (1) be a true enough model and will 
be used to evaluate the designed controllers. 

It is understood that the triple matrices ( M e ,C e ,K e ) in (1) do not have the usual structural properties, 
that is, they are neither symmetric nor sign definite. It is apparent from (1) that the aircraft rigid body 
dynamics and the aeroelastic modes are coupled. When the couplings are neglected, the performance of flex- 
ible modes and aircraft open- loop behavior can be determined by examining the triple matrices (M e , C e , K e ) 
and the eigenvalues of A a . We regard this as ’’‘uncoupled’” condition, which shows some underline system 
characteristics. The detail derivations and physical interpretations of (1) can be found, for instance, in Ref. 
6 . 

For better illustration, we rewrite (1) in a state-space representation form as follows, 

Xp = Ap Xp T Bp us , (2) 

where x p = (x e , x e , x a ) and us = ( S e , 6f, S s ), and 



1 

0 

I 

1 

O 


l 

O 

O 

O 

l 

A.p — 

0^ 

—M~ 1 C e 

—M~ 1 C a 

II 

of 

0 M-'D f M-'D S 


CO 

1 

Ad 

Aa 


B e Bf B s 


We assume that the pair (A p , B p ) is controllable. In this paper, we will consider modal reduction process for 
flexible wings, hence it makes little practical sense if the residual modes are unstable. Therefore, we assume 
A p is Hurwitz and has distinctive eigenvalues. 


II. A. Wind Turbulence Model 

We consider the coupled aeroelastic aircraft model given in (2) is subject to a turbulent wind gust load. In 
this analysis, we assume that the wind gust model is described by 

I A,, — A w x w T B w Wg 
Dw C w x w T D w Wg 

where x w denotes the states, w g the random gust wind which is assumed to be zero-mean white noise with 
intensity W, and y w the total random gust load applied to both rigid body aircraft and aeroelastic wings. 
We assume that A w is Hurwitz. If Dryden’s longitudinal wind turbulence model is used, then x w will be 
consisting of vertical airspeed, vertical acceleration, and pitch rate. Adding the gust load y w to (2) to obtain 

x p = ApX p + B p us + y w . (4) 
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II. B. Actuator Dynamics 


The actuator dynamics that produce control surface deflections us are assumed to be derived from the 
following linear systems, 

( xg = A s xs + B s u ^ 

[ Us = CsXs 

where Xs is the states of actuator dynamics and u the control command. It is assumed that the pair (A^, Bs) 
is controllable and As Hurwitz. For a 2nd order actuator, xs will be consisting of control surface deflection 
angle and angular rate. It is important to note that, in practice, xs is bounded. 


II. C. Open-loop System Representation 

After substituting (3) and (5) into (4), we obtain the overall open- loop state-space representation as follows, 


j] = A v r] + B v u + D v w g (6) 

where ry = ( x p , x w , xs) and 
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It should be noted that the turbulence state x w is not controllable. In this paper, the control/performance 
outputs are denoted by y and the sensor measured outputs by z, and they are described by 


\ y = c ,v (7) 

[ 2 = E v r] + v 

where C indicates the control outputs of interest and E v the locations of measurement sensors, and v 
denotes sensor noise which is assumed to be zero mean white noise with intensity V. We can partition the 
performance output y into a series of block outputs as 


V i 
Vi 

Um 


Vi = {C v )iV . 


( 8 ) 


The reason for grouping control output y in (8) is that for systems of large dimensions, the control objectives 
and constraints are often prescribed with respect to a collection of system states. For instance, we may 
choose y\ = xs to represent the constraints on actuator position and/or rate limits, and (t/2, • • • , y m ) the 
constraints on deflections and/or rates for aeroelastic wings. 

Combining (6) and (7) yields the following full-order open- loop system description, 


£ : 


y 

y 


z 


A v y + B v u + D v w g 
C v y 

ErjT] + V 


(9) 


Note that the dimension of 77 can be very large. Henceforth, (9) or its equivalent form is called full-order 
evaluation model, while the design model, to be introduced later, represents the reduced-order model from 
which a stabilizing controller is designed. 


II. D. Control Design Objective 

The aircraft flight control system development cycle often requires optimization of control law subject to a 
chosen set of handling qualities specifications.' In this study, we let yspec to denote such set, which captures 
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aircraft rigid body performance, and let Si to denote the Level 1 requirements for yspec- For controlling 
the vibrational motion of aeroelastic wings, we let Yj > 0 denote the desired output covariance constraint 
matrix for yj, where j = 1, 2, • ■ • , m. In this paper, the control design objective is to find a stable, strictly 
proper, dynamic output feedback controller of the form 


X c — A c £ c d" L c z 
U = K C X C 


, x c eR n 


(10) 


that stabilizes the full-order evaluation model £ described in Eq. (9), while minimizing the weighted control 
effort 

U = lim £{u*(t)-Ru(t)} , (11) 

t — yoo 

subject to 

1. yspec € Si or all specifications are in Level 1 regions, and 

2. output covariance matrix Yj for aeroelastic wings is constrained by Yj, i.e. 

lim £{yj(T)yj(T)} = Yj < Yj , j = 1, 2, • • • , m . (12) 

t — J 


Note that L c and K c are control gain matrices, £ denotes an expectation operator, and R is a sign definite 
symmetric weighting matrix. In general, the order of evaluation model can be very large, hence it is neither 
practical nor necessary to implement a full-order controller for (9). However, one may choose to design a 
full-order controller for initial assessment of reasonable or achievable output covariance constraint matrices 
of the closed-loop system, and use these to set Yj for design iterations. In this paper, we will design a 
controller of the form (10) for a reduced-order design model. 



Figure 1. Closed-loop block diagram for aeroelastic aircraft model 


III. Modal Cost Analysis 

Beacuse of practical limitations on control bandwidth and to reduce computational burden in actual 
implementation, model reduction is often an essential part of controller design for physical systems of large 
dimension. This paper adopts the notion of modal cost analysis (MCA) 5 and utilizes it to attain a reduced- 
order design model. The concept of MCA was derived from the component cost analysis (CCA) , which was 
first introduced in Ref. 4 and applied to the control of large space flexible structures. 4 The idea of CCA is 
to examine the contribution, and hence significance, of each state component to the mission objectives in a 
control system. A metric of component contribution can be calculated in terms of output covariance, from 
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which contribution of each component was studied and ranked from high to low. The approach was used to 
derive a reduced-order model by neglecting less significant components. The MCA is an extension of CCA to 
the standard nratrix-second-order systems, where mass, damping, and stiffness matrices are sign definite and 
symmetric. A modal transformation can be applied so that in modal coordinate each mode is decoupled. In 
this case, contribution of each mode to output covariance is examined directly and a reduced-order model can 
be determined by retaining only dominating modes. Applications of MCA to model reduction and control 
system design was studied in Refs. 8 and 10. In this paper, because of the coupling between aerodynamics 
and aircraft dynamics, the total mass, damping, and stiffness matrices are neither sign definite nor symmetric, 
as indicated in (1). However, the MCA method can still be applied with minor modification. 

Recall £ described in (9), 

! fi = A v ij + B v u + D v w g 

y = c vV 

z = E v ri + v 

It should be noted that the modal cost analysis is performed only to the aeroelastic modes. Since A p is 
Hurwitz and has distinctive eigenvalues, there exists a nonsingular transformation matrix T such that, with 
rj = Tfj , £ can be transformed into a block decomposed form as 


where f] is partitioned into 


t 


f) = At) + Bu + Dw g 
y = Cfj 
z = Efj + v 


m 

m 


v = 


y2N 


(13) 


and fji € R 2 represents the modal states for Ah aeroelastic mode; i = 1, 2, • • • , 2N. Accordingly, the system 
matrices can be partitioned as follows, 


A = T 1 A V T = diag{ Ji, J 2 , • • • , J- 2 N, • • ■ } I A 


di ~h 

bi a, ; 


B = T~ 1 B rt = 


B x 

B 2 

B 2 n 


,C=Cr,T = 


Cl C 2 


C 2 N 


D = T~ 1 D r> = 


D x 

B>2 

D 2 n 


E = EgT = 


E 1 E 2 ■ ■ ■ E 2 n 


(14) 


Note that the matrix A denotes a real block diagonal form for a pair of complex eigenvalues a; ± jbi. 
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III. A. Computing Modal Cost 

The MCA considered in this paper is to study modal contribution of each aeroelastic mode, denoted as zy, 
to a weighted output covariance cost function V given by 

2N 

V = lirn 8 {V (t)Q?/(t)} = V] , (15) 

i=l 

when subject to random input w g . In (15), Q is a block diagonal symmetric and non-negative matrix, and 
r'o denotes modal cost contribution from other modes such as aircraft rigid body dynamics and actuator 
dynamics, etc. In what follows, we calculate the open- loop modal cost jy. 4 

Recall (13). Since A v is Hurwitz and the pair (A V ,D V ) is stabilizable, the following Lyapunov equation 

PA t + AP + DWb t = 0 (16) 

renders a unique positive semi-definite symmetric solution P, where we recall W is the covariance of w g . 
Next, the modal cost i y can be obtained by 4 

u^tr^QCPU, i = 1,2,--- , 2N , (17) 

where tr(X) denotes the matrix trace of X and [ (PQCP ]„ is a 2 x 2 output covariance matrix corresponding 
to ith aeroelastic mode. Now, given iy we then can rank them based on their modal cost as 

N > |^| > •••Kb" > \"2n\, ( 18 ) 

where iq being the most critical mode and v 2 n least critical mode. Note that zy can be negative, which 
indicates that this particular mode is in fact helping to reduce the total cost V, 4 however the total cost is 
non-negative. It should cause no confusion that in (18) zy means the ith rank modal cost, not necessarily 
the modal cost of ith aeroelastic mode. It is important to emphasize that the choice of weighting matrix Q 
directly affects the modal cost. 

Remark: When the open-loop system contains unstable mode(s), the proposed MCA model reduction 
process can still be applied with modification. Since the covariance cost for unstable modes is infinity, these 
modes should by default be retained. The MCA is then performed for stable modes, following the same 
process as presented. 

III.B. Reduced-order Design Model 

Based on (18) we may choose to keep the significant n aeroelastic modes, where n <C 2 N. Furthermore, n is 
the number of modes that can be accommodated in the control system synthesis, given its bandwidth. If n 
is chosen, then we may re-arrange the modal states f) in (13) and decompose it into ( x r ,Xt ), where x r and 
Xt denote the retained and truncated modal states, respectively. Therefore, (13) can be rewritten as follows, 



Recall that the retained state x r consists of n aeroelastic modes, aircraft rigid body dynamics, actuator 
dynamics, and turbulence model. Now, we can present the following reduced-order design model to which a 
stabilizing controller will be designed. 

! x r = A r x r + B r u + D r w g 

y r = C r x r (20) 

z r = E r x r + v 
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In £ r we have neglected the influence of observation spillovers induced by the truncated high-order aeroelastic 
modes. Since the modal cost contributions of these modes to the output covariance are less significant, it is 
expected that the reduced-order design model £ r captures the salient feathers of the full-order model, and 
hence is used to build a dynamic output feedback controller. Ultimately, this controller will be evaluated 
with full-order evaluation model (9). 

To summarize the open-loop model reduction process, we first need to perform coordinate transformation 
to (9) to block diagonalize the system matrices. Then, solve the Lyapunov equation (16) for P and calculate 
the modal cost i>i, and finally rank v,; and determine the number of aeroelastic modes to be retained. 


IV. Integrated Flight Control and Output Covariance Constraint Problem 

Recall the dynamic output feedback controller given in (10), 

J ■ 1 , 1 — A c x c T C c z 
\ U = K c x c 

Here, we choose A c = A r — L c E r + B r K c and x c G R Ur . Though the order of (10) is chosen to be the 
same as that of reduced-order design model, it utilizes the full measurement information, including the 
truncated modes. When interconnecting (19) with (10), we obatin the feedback-controlled full-order system 
representation as 

| i = AX + VW 

| y = cx 


(21) 


where 
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Note that in A the upper left 2x2 block matrix is precisely the closed-loop matrix when closing (10) with 
the design model £ r , and the off-diagonal terms are due to observation and control spillovers of truncated 
modal states. Therefore, the design of control gains K c and L c is first performed for the design model £ r . 

IV. A. Output Covariance Control for Design Model 

In this section, we consider the following reduced-order closed-loop system representation which is extracted 
from (21), 

f X, = A r X r +V r W 

\ y r = c r x r 

where 

r 

Xf' x c 

X c 


(22) 


V. = 


, y r = 


Vr 
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and 


A r = 


./4'p J-J Q 

0 

, v r = 
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L c 
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K c _ 


The output covariance constraint (OCC) problem 9 for (22) is stated as follows: 
For £ r , find a dynamic output feedback controller of the form 


x c — ( A r — L c E r T H r A c )x c T L c z r 

u = K c x c 


(23) 
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so that it minimizes the control cost 


lim £ {u* (t) Ru(t)} ; R > 0 , 


subject to the output covariance constraint matrices Y,. i.e. 

{C r ) i P{C r )\=Y i <Y i -i = 1,2,--- ,m, 
where P > 0 solves the following Lyapunov equation 

A r P + PA* + D r WD* r = 0 , 
and W = diag (W, V) denotes the covariance for W. 

The following lemma provides the solution to the OCC problem for (22). 

Lemma l. 9 Consider the closed-loop system defined in (22). Suppose ( L C ,K C ) is the optimal solution to 
the OCC problem. Then, there exist a positive definite symmetric matrix R and a block diagonal matrix Q 
of the form 

Q = diag {Q i,Q 2 ,- ■■ ,Qm} > 0, (24) 

such that the control gain matrices K c and L c are given by 

K c = —R~ 1 B t r X , L c = YElV - 1 , (25) 

where X > 0 and Y > 0 are solutions to the control and filtering Riccati equations, 

0 = XA r + A\X - XB r R~ 1 B t r X + C* r QC r 
0 = A, Y + YA t r - YE*V~ 1 E r Y + D r WD* r 

Furthermore, the chosen weighting matrices and control gains satisfy the following output covariance con- 
straints, 

0 =(Yi- Yi ) Qi 

Yi = (C r )i(Y + X c ){C r )\ ; * = 1, 2, - - ■ ,m, 

where X c > 0 solves 

0 = (A r + B r K c )X c + X c (A r + B r K c y + L c VLl . 

Therefore, if K c and L c are chosen according to (25), then both (A r + B r K c ) and (A r — L c E r ) are Hurwitz, 
so are A r and A. since A t is Hurwitz. In Ref. 9, an algorithm for selecting weighting matrices was presented 
and an optimal solution to the OCC control problem obtained, as stated in the Lemma. It is important 
to note that the OCC problem described above is for the reduced-order design model. We will incorporate 
this process into the OCC control problem for full-order evaluation model when we develop an integrated 
modeling and control approach in Section V. 

IV. B. Optimal Flight Control Design in CONDUIT 

CONDUIT (CONtrol Designer’s Unified InTerface) is a state-of-the-art multidisciplinary computational fa- 
cility for aircraft flight control design, evaluation, and integration. It incorporates aircraft dynamic models, 
handling qualities evaluation, closed-loop system performance, and multi-objective optimization into a MAT- 
LAB/SIMULINK based interactive environment. For detailed description of CONDUIT and its applications 
to flight control development projects, one may refer to Refs. 11-14. Here we provide only a brief overview 
of CONDUIT. 

To setup a problem in CONDUIT, we first need to define the aircraft dynamics, actuator dynamics, and 
architecture of the control law in SIMULINK. Then, select a set of appropriate design specifications for 
performance evaluation from the provided specification libraries, or create user-defined specifications using 
CONDUIT’S SpecMaker tool. In each design specification there are three levels of compliance, namely, Level 
1, Level 2, and Level 3. 12 A sample design specification for stability gain/phase margin (StbMgGl) is shown 
in Figure 2, where Level 1 region, colored in blue, is the most desirable performance region indicating gain 
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StbMgGI : Gain/Phase Margins 
(rigid-body fteq. range) 
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Figure 2. Sample Design Specification: Stability Gain/Phase Margins 


margin above 6 dB and phase margin above 45°. Level 2, colored in magenta, means the performance is 
adequate, and Level 3, colored in red, indicates an inadequate performance region. Note that Level 1/2/3 
boundaries can be adjusted to suit user’s specifications. 

Next, prioritize each design specification by designating it as hard constraint, soft constraint, or summed- 
objective. 11 Finally, CONDUIT utilizes CONSOL-OPTCAD as optimization engine and conducts multi- 
objective optimization in three distinct phases. In Phase 1, the control law is tuned so that all hard 
constraints reach Level 1 region. Phase 2 will work on the soft constraints and ensure they reside in Level 
1 region, whereas in Phase 3, CONDUIT will optimize the control law to the selected summed-objective 
performance criteria, therefore to attain the best control design from the family of feasible solutions. 

In this paper, the design specifications include not only for aircraft rigid body dynamics, but for aeroelastic 
wings. The vibrational behavior of the wings, which is captured in yj , will be constrained as output covariance 
Yj and needs to be incorporated in CONDUIT as part of handling qualities specifications. Therefore, we 
need to create a set of OCC specifications for full-order evaluation model (9). In what follows, we assume 
yj be a scalar output taken at some location along the flexible wing, and Yj a desired output covariance. 
Let Yj > Yj be an acceptable output covariance for yj. Then, using SpecMaker tool in CONDUIT, we can 
create an OCC design specification, AseOcGl, shown in Figure 3. 

Therefore, in CONDUIT the complete handling qualities and performance specifications will consist of 
those for rigid body dynamics, yspecj and aeroelasticity. CONDUIT optimization process will then tune the 
weighting matrices, and hence the control gains ( K C ,L C ), until j/spec G Si and Yj < Yj. 

V. Iterative Modeling and Control Design Process 

We present in this section an iterative algorithm that involves model reduction, OCC controller design, 
and optimal flight control design, subject to aircraft handling qualities and performance specifications. The 
design process searches for the controller with the best performance by tuning the design model until the 
iteration converges, as illustrated in Figure 4. The model reduction utilizes MCA, and the design of flight 
controller for rigid body dynamics and aeroelastic wings follows the procedure presented earlier. This is 
essentially a closed-loop model reduction process, since the choice of weighting matrix Q directly affects the 
modal cost, and hence model reduction process. On the other hand, matrix Q also affects the computation 
for output covariance, in that if a particular output covariance specification Yj is not satisfied then the 
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AseOcGI : Generic OCC spec 



Figure 3. OCC Design Specification 


corresponding Qj will be adjusted for next iteration. Therefore, by combining model reduction and OCC 
control problem, we are able to find an appropriate weighting matrix Q that reflects the importance of 
each output and that the reduced-order design model retains the dominant dynamic characteristics which is 
important to attain the required flight control performance. 

CONDUIT’S external scripting capability 14 is quite suited for performing the iterative control design 
process depicted in Figure 4. The external scripting file is a MATLAB rn-file which allows us to incorporate 
the OCC control problem in CONDUIT and executed with optimization process. 

VI. Applications to Generic Transport Model 

Here, we present an application of proposed integrated OCC algorithm to longitudinal aeroelastic GTM 
as an illustraive example. GTM is configured with 6 trailing edge flaps and 6 leading edge slats, in addition 
to the elevator, see Figure 5. The objective is to design an optimal reduced-order output feedback controller 
for ” ‘full-order’” GTM, so that the closed-loop system satisfies the desired handling qualities/performance 
requirements, and the vibartional motion of the wings is suppressed. 

In this paper, GTM at cruise condition at Mach 0.8 and altitude at 30, 000 ft with 50% fuel remaining 
is considered. Furthermore, we let N = 10, that is, we account for 10 bending and 10 torsional modes. It is 
assumed that both longitudinal rigid body states and aeroelastic wing states are available for measurement. 
For aeroelastic measurements, we take the outputs at 10 equally spaced data points along the wing span; 
and the last point being at the wing tip. We measure both displacement and rate at these 10 points. The 5 
control outputs for aeroelastic wings are assumed to coincide with the last 5 measurement points. For this 
analysis, we assume that GTM is subject to a random gust turbulence of ±20 ft/s and random measurement 
noise of 0.001ft 2 (rad 2 ). 

VI. A. Open-loop Analysis for GTM 

We choose Q to be an identity matrix and follow the open-loop MCA procedure described in Section III. 
Table 1 contains the list of modal costs, along with their modal frequency and damping. 

The open-loop MCA shown in Table 1 indicates that beyond 4th bending and 2nd torsional mode the modal 
cost becomes less significant. Higher mode demands larger control bandwidth, therefore, given the frequency 
range of actuator dynamics and results of MCA, we can decide specific aeroelastic modes to retain in the 
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Figure 4. Integrated Modeling and Control Algorithm 


Table 1. Open-loop Modal Cost Analysis for GTM 


Mode ID 

Modal Cost (10 6 ) 

Frequency (rad/sec) 

Damping 

1st torsion 

0.281 

8.47 

0.011 

1st bending 

0.017 

8.71 

0.289 

2nd bending 

0.039 

14.52 

0.210 

2nd torsion 

1.326 

15.70 

0.007 

3rd bending 

0.0027 

24.49 

0.077 

4th bending 

0.0024 

31.86 

0.067 
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(Re-print from IPP Report by N. Nguyen) 


Figure 5. Generic Transport Model 
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design model. We may start with retaining 6 modes; namely, 1st to 4th bending mode and 1st to 2nd 
torsional mode, and design a reduced-order dynamic output feedback controller given in (10) by following 
the procedure presented in Section IV. A. 

Before proceeding, we need to access the level of output covariance of the open-loop system. This will help 
us gauge and setup the achievable output covariance level for closed-loop system for optimization process. 
The output covariances for open-loop system are shown in Table 2. 

Table 2. Open-loop Output Covariance 


Data Point 

Bending Deflection x 10 2 (ft 2 ) 

Torsional Deflection x 10 6 (rad 2 ) 

Pt. 1 

0.52 

1.71 

Pt. 2 

1.03 

2.79 

Pt. 3 

1.71 

4.07 

Pt. 4 

2.52 

5.38 

Pt. 5 

4.17 

6.79 


VI. B. P rob elm Setup in CONDUIT 

The emphasis in this section is to define an iterative control optimization problem to be solved in CON- 
DUIT. We have chosen six design specifications for analysis; they are: Eignevalues, Stability Gain/Phase 
Margins, Quickness, Crossover Frequency, Actuator Saturation, and Normalized Output Covariance. The 
eigenvalues specification verifies the closed-loop system is stable, the stability margin specification ensures 
that the satisfactory gain and phase margins are attained for the broken-loop response, the crossover fre- 
quency specification is to ensure minimal overdesign, pitch axis quickness low bandwidth specification is 
to ensure the ratio of peak attitude rate to change in attitude be exceeding specified limits, the actuator 
saturation specification is used to limit the actuator saturation in position and rate over a specified period, 
and the output covariance specification is to ensure the vibrational motion of the wings is contained. These 
specifications are listed in Table 3, along with their constraint type and level I specifications. 

Table 3. Design Specifications Used in CONDUIT 


Requirement 

Specification 

Source 

Constraint Type 

Level I 

Eigenvalues 

EigLcGl 

Ames Research Center 

Hard 

Re(Ai) < 0 

Stability Margins 

StbMgGl 

MIL-F-9490D 

Hard 

at least 6dB/45° 

Quickness 

QikAtGl 

Ames Research Center 

Soft 

As shown 

Crossover Frequency 

CrsLnGl 

Ames Research Center 

Summed Objective 

< 2 rad/s 

Actuator Saturation 

SatAcGl 

Ames Research Center 

Soft 

1 sec 

Output Covariance 

AseOcGl 

Ames Research Center 

Hard/ Soft 

< 10% x OC open 


Figure 6 shows the handling qualities (HQs) window in CONDUIT. It should be emphasized that the 
objective of OCC problem for aeroelastic wings is to achieve an order of magnitude less than the output 
covariance of the open- loop case, as indicated in Table 3. In this regard, Level I/Level II boundary of 
AseOcGl is normalized to 1. Note that there is a built in scale factor within specification that we could use 
to re-scale the level I region, if needed. 

The tuning design parameters are contained in the ’Q’ matrix. Let dpp_sp to indicate the tunable 
weighting parameter for longitudinal short period mode, dpp_bend for bending deflections, and dpp_tors for 
torsional deflections. It should be noted that, since the aircraft’s rigid body dynamics is coupled with aeroel- 
stic modes, and in addition the bending and torsional motions are coupled, the three weighting parameters 
can not be tuned independently. 
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VI.C. CONDUIT Optimization 


In this section, an optimal dynamic output feedback controller is designed for the reduced-order design 
model, which consists of rigid body short period mode and 6 aeroelastic modes listed in Table 1. The 
controller is developed by solving the two Riccati equations given in Lemma 1 with varying weighting matrix 
Q. Subsequently, this reduced-order controller is applied to the full-order system S in (9). The iterative 
optimal selection of weighting matrix is performed in CONDUIT subject to the list of design specifications 
given in Table 3. Figure 6 shows the evaluation of HQs with initial seletion of weighting parameters, i.e. 
dpp_sp = dpp_bend = dpp.tors = 0.0001. These parameters reflect a small control effort by closing the 
feedback loop for initial check up. The results show that the closed-loop system is stable and satisfies overall 
performance requirements, but the output covariances in bending and torsional displacements fail to meet 
level I requirement. Figure 7 shows the HQs of optimal solutions after 19 iterations. It shows that all design 
specifications reach level I regions, and the optimal weighting parameters are clpp_sp = 0.2127, dpp_bend 
= 0.0021, dpp.tors = 17.6356. Figures 8 and 9 show the comparisons between open-loop and closed-loop 
bending and torsional deflections at Pt.l and Pt.5, respectively, subject to random gust turbulence. 

We have also designed optimal output feedback controllers for varying number of retaining aeroelastic 
modes; namely, 5 modes (B1,B2,B3,T1,T2), 4 modes (B1,B2,T1,T2), and 3 modes (B1,B2,T2). In all these 
cases, the closed-loop full order system is stable and satisfies relevant performance requirements, except 
the output covariances, especially, the torsional deflections at outer wing (Pt.5). Figures 10 and 11 show 
the closed-loop output covariance comparisons with various design models. As expected, the closed-loop 
performance degrades as the order of controller reduces. 

VI. D. Flutter Analysis 

Figure 12 shows the root locus of the open-loop aeroelastic model as aircraft speed varies from Mach 0.6 
to 0.9, and at around Mach 0.82 the system becomes unstable (B1,T1). The designed optimal reduced- 
order controller is used for flutter analysis. Figure 13 shows the root locus of the closed-loop system with 
varying air speed from Mach 0.8 to 0.92, and Figures 14 and 15 show the output covariance performance 
comparisons. The closed-loop system is the interconnection of the open-loop aeroelastic model with the 
proposed reduced-order dynamic output feedback controller based on 6 aeroelastic modes. Note that the 
proposed dynamic controller is designed at Mach 0.8, the results in these figures indicate that the designed 
controller is fairly robust against modeling uncertainties, though the performance degrades somewhat. As 
expected, the closed-loop system stability breaks down when air speed continues to increase, in this case, 
beyond Mach 0.88. 


VII. Conclusion 

In this paper, we have presented an integrated model reduction and control system design process for 
aircraft with flexible wings. Model reduction is conducted using modal cost analysis for each aeroelastic mode, 
and the modes that show significant contribution to the output covariance are retained. Since the notion of 
output covariance control is new to the aircraft aeroelastic study, in this paper we have used 10% of the open- 
loop output covariances as the closed-loop design requirement for aeroelastic modes. The optimal dynamic 
output feedback controller is designed for the reduced-order design model, by following a procedure similiar 
to L 2 / -Hoo multiobjective control design process. The process is to find the proper weighting parameters 
iteratively until the design objectives are met. The output covariance control algorithm is implemented in 
CONDUIT for full-order evaluation model. The analysis shows that the proposed controller is able to achieve 
level I design specifications for both rigid-body aircraft performance and flexible wings. It is also shown that 
the proposed controller is robust in the presence of modeling uncertainties. 
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Figure 6. Handling Qualities at baseline gains 
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Figure 7. Handling Qualities at optimal gains 
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Figure 8. Bending displacement comparison between open and closed-loop at Pt. 1 and Pt. 5 
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Figure 9. Torsional displacement comparison between open and closed-loop at Pt. 1 and Pt. 5 
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Output Covariance (ft) 



Figure 10. Output covariance (bending) comparison for 6-mode, 5-mode, 4-mode, and 3-mode 



Figure 11. Output covariance (torsion) comparison for 6- mode, 5- mode, 4-mode, and 3- mode 
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Figure 12. Root locus as air speed varies from Mach 0.6 to 0.9 
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Figure 13. Closed-loop root locus as air speed varies from Mach 0.8 to 0.92 
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Figure 14. Closed-loop output covariance (bending) comparisons 



Figure 15. Closed-loop output covariance (torsion) comparisons 
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